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ABSTRACT 

We establish the existence of concave solutions of Berman's equation 
which describes the laminar flow in channels with injection through porous 
walls. It is found that the (unique) concave solutions exist for all 
injection Reynolds number R < 0. 
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I. INTRODUCTION 


Studies of the transpiration- or sweat-cooling of heated surfaces such as 
turbine blades, rocket walls, or wing surfaces in high-speed flight by the 
diffusion of fluid through porous-metal combustion-chamber liners have been 
made within the framework of the boundary-layer flow along a porous plate with 
fluid injection. Laminar flow in a two-dimensional rectangular channel with 
porous walls was studied by Berman in 1953. He showed that the corresponding 
Navier— Stokes equations can be reduced to a nonlinear third order ordinary 
equation with two point boundary conditions and Reynolds number R based on 
injection-velocity. Berman also gave perturbation results for extremely 
small R. In later years, various asymptotic results and numerical results 
(e.g., Sellars (1955), Yuan and Finkelstein (1956), Terrill (1964), Raithby 
(1971), and Robinson (1976)) were given. Asymptotic results have of necessity 
been purely formal, since the existence of a solution has not been proved nor 
the necessary estimates obtained. 

Skalak and Wang (1978) classified all possible solutions of Berman's 
equation, both for injection and suction. It was shown that there is only one 
possible solution for the injection problem. In this note we establish the 
existence of this solution. In Section 2, we derive the Berman's equation. 
In Section 3, we give some preliminary results and several important a-priori 
estimates of the solutions. In particular, we show that the maximum axial 
velocity in the channel is insensitive to the Reynolds number R. In Section 
4, we use the Learay-Schauder Fixed Point Theorem to establish the existence 
of solutions. In Section 5, we give some numerical results which confirm our 
a-priori estimates and agree with the asymptotic results obtained by Yuan and 
Finkelstein (1955) as R ->• -°°. 
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II* DERIVATION OF BERMAN'S EQUATION 

Consider the steady, incompressible, laminar flow along a two-dimensional 
channel with porous walls through which fluid is injected with uniform speed 
v w . Take x and y to be coordinate axes parallel and perpendicular to the 
channel walls and assume u and v are the velocity components in the x 
and y directions respectively. Letting the channel width be 2h and 
introducing the dimensionless variable 


n 



reduce the Navier-Stokes equations to 


(i) 


( 2 ) 


u 


2 2 

3u u 3u _ 1 3p ,3 u 1 3 Uv 

3x h 3n ~ F 3x \2 72 TT' 

3x h 3n 


3v + v 9v_ 
3x h 3n” 


l_ 3p 
hp 3y 


+ v 


2 2 
,3 v 1 3 v. 

' 2 + 2 2 ^ 
3x 3n 


where p, p, and v are the density, pressure, and kinematic viscosity 
of the fluid respectively. 

The continuity equation is 


(3) 


3u _1 3v 
3x h 3y 


= 0 . 


With the walls of the channel at y = ±h; i.e., q = ±1 , the boundary 


conditions are 
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u(x, ±1) = 0; v(x,0) = 0, 

(A) 

v(x, ±1) = constant = v <0; (x,0) = 0. 

’ w 9y 

For a two-dimensional incompressible flow, a stream function ij> exists 
such that 


(5) u(x ’ n) = If 

(6) v(x, n) = - 


and the continuity equation, (3), is satisfied. 

Since the flow is symmetrical about a plane midway between the walls, the 
solution will be investigated over half the channel, i.e., from the midplane 
to one wall. 

A suitable stream function is 

(7) iKx,A) = [hu(0) - v x] f(n) 

w 

where u(0) is the entrance velocity at x = 0* Thus the velocity 
components u and v are given by 

(8) u(x, q) = h ^huCO) - v x] f'(n), 

w 


v(x,n) = v f(n). 

W 


( 9 ) 
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Substituting (8), (9) into (1) and (2), we have 


( 10 ) 



9x3ri 


0 


( 11 ) 


+ R(f ' 2 - ff~) = K 


where R = v h/v <0 is taken to be the injection Reynolds number of the 
w 

flow and K is an integration constant, to be determined later as a function 
of R. The boundary conditions (4) become 

(12) f(0) = 0, f"(0) = 0, 

(13) f'(l) = 0, f(l) = 1. 


Thus, the nonlinear two-point boundary value problem we are considering 
consists of (11), (12), and (13) with R < 0. By a solution of (11), (12), 

3 

and (13) we mean a function f € C [0, 1] and a constant K which satisfy 
(11), (12), and (13). 

Remark : It is easy to see that once we have a [0, 1] solution f, then 

f e C k [0, 1] , k 2 0. 


III. PRELIMINARY RESULTS 

If we differentiate (11), we get 

f (iv) 


(14) 


+ R(f'f" - ff'~) = 0 



Lemma 15: JLf f Is a solution of (11), (12), and (13), then f 

is increasing on [0, 1] . 

Proof ; It is sufficient to show that f^ iV \n) 2 0 on [0* 1]. 

Differentiating (14) again, we get 

(16) f* v * = - Rf" 2 + Rff( iv ). 

From boundary conditions (12), f^ iv ^(0) = 0. (16) shows that at any point 

c where f (iv) (c) = 0, f (v) (c) = -Rf" 2 (c) >0. Hence, f (iv) ( n ) > 0 

Vn e [0, 1]. 

Remark . We have shown also that f'f" - ft"' > 0 Vn [0, 1], 


Lemma 17: Let f'(0) = a, f'"(0) = $. Then for R < 0, 

(2 1 3 < "3; -3 < f"(l) < - | 

(ii) 0 < f'(n) < a; 3 < f — (n) < 0; f — CD < f"(x) < 0, on [0, l]. 


Proof_: Let G 4 (n,t) be the Green's function of 


(18) 


f (iv) = 0; 


and 


(19) 


v(0) = v" ( 0 ) = v'(l) = v( 1 ) = 0, 
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which is given by 

( nCt-i) 2 [(3 - n 2 )t - 2n 2 ], 0 < n < t, 

t(n - 1) 2 [(3 - t 2 )n - 2t 2 ], t < n < l. 

Then the solution of (14) can be written as 


(21) 

f(n) = R / 1 G(n , t)(ff'~ - f'f")dt + f Q (n) 

where 


(22) 

f Q (n) = f (3 - n 2 ). 

From (20), 

we can write 

(23) 

a = I + T R /^(i - t) 2 ( ff '" - f ' f ") dt 
Z 0 

(24) 

B = -3 - | / 1 (2 + t)(l - t) 2 (ff'~ - f'f'")dt 
0 

(25) 

f"(l) = -3 + y J^tU - t 2 )(ff"'- f'f")dt. 
0 


By the remark following (16), we conclude that 
(26) o > | , B < -3, f~(l) > -3. 
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From Lemma 15 and equation (11), we have 

f (0) = S = K - Ret 2 < f (1) = K + Rf"(l). 


Thus , 

< 27 ) f"(l) < -a 2 < - | < 0. 

Combining (26) and (27), we get (i). Now we will prove that f~(n) is 

rnonotonically decreasing on [0, 1], 

From (20) and (21), it is easy to show that 

(28) f'~(n) <0 Vn e [0, l]. 

Indeed, should f (c) = 0 at Vc [0, 1], Assume c is the first such 
point; then 

f (iv) (c) = (-R)f 'f "(c) < o. 


which contradicts the fact that f( lv ) > q on [0, 1], Hence, 

f (b) K 0 and both f " and f ' are rnonotonically decreasing on [0, 1]. 
f ' looks like 
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Now let us get a rough lower bound for g . By the integral 
representation (24), since (2 + t) ( 1 - t 2 ) _< 3 for t £ [0, 1], we have 

(29) 6 > -3 - | R/^ff'" - f'f"")dt = -3 - 2 R(f"(l) + a 2 ). 

0 L 

We can get other lower bounds on g, e.g., by integrating (11) once and use 
(12), (13) and the fact that 

(30) K = 0 + Ra 2 
we get 

(31) 0 > -2a + R(a 2 - 1). 

We prefer (29) to (31), because from the asymptotic expansions of 
f (n) as R + -°° given by Yuan and Finkelstein (1956), we have 

f~(l) + a 2 = 0(4-) as R -V 

K 

then by (29) and the fact that g is monotone decreasing with respect to 

|r|» one can expect that g is bounded by a fixed number (which turns out 
3 

to be - — , see Section 5). 

IV. EXISTENCE THEOREM 

Let {E = u€C 2 [0, 1], u(0) = u~(0) = u'(l) = 0} . 

space if E is given the usual sup-norm of C^[0, 1], i.e., 


E is a Banach 
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(32) Hull = Sup |u(n)| + sup |u'(n)| + sup |u"(n)| 

n e [o, 1] n e [0, l] n € to, l] 


Our system (11), (12), and (13) is equivalent to 


(33) 


f(n) = h(n) + R/ 1 [J(n, t) - h(n)J(l, t)](f' 2 - ff")dt, 


and 


(34) 


K = -3 + SRj^JO, t)(f' 2 - ff")dt, 
0 


where J(n, t) is the Green's function of -v"' = 0, v(0) = v"(0) = 
v'(l) = 0 and is given by 


(35) J(n , t) = 


(1 - t)n 


i 1 f 2 2 v 

- 7 (n + t ) + n 


0 < n < t, 


t < n < l; 


while 


h(n) = 3/ 1 J(n, s)ds. 

0 

Note that (33) is independent of K. Any solution f of (33) in E will 
automatically satisfy f(l) = 1, since h(l) = 1, 

We will use the following version of the Learay-Schauder Fixed Point 
Theorem to prove the existence of solutions to (33): 
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Theorem 36 (Learay and Schauder (1934)): Assume that the operator 

^( u » X) (ueE, 0 <_ X 1) are completely continuous and A(u, 1) = Au. 

Assume that all solutions of the equation u = A(u, X) (0 < X < 1) satisfy 

the common a priori bound null < r Q . If the vector filed u - A(u, 0) has 
nonzero Learay-Schauder degree on spheres {u|llull = r} of larger radii 
r > r Q> then the equation (I - A) u = 0 has at least one solution. 

We shall define our operator A(u, X) in the following way. Let 

(37) Lu = u - h. 

Then L _1 exists and equation (33) is equivalent to 

(38) u = A(u, X) = L _1 [XR/ 1 [J(n, t) - h(q)J(l, t)](u' 2 - uu")dt] . 

0 

To show that A(u, X) is completely continuous on E, since L -1 is linear 
and bounded, it is sufficient to show that the integral operator on the right 

of (38) is compact from E E for all X e [0, 1]. This is easily shown 

by noting that the operator 

(39) F(u) = u' 2 - uu~ 

is a continuous bounded mapping from E •> C°[0, 1], i.e., F maps bounded 
sets of E into bounded sets of C®[0, 1], and that the operator 

Bv = / 1 [J(n, t) - h(n)j(l, t) ] v dt 
0 


( 40 ) 
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is compact from C^[0, 1] + E. A simple application of Arzela-Ascoli 
Theorem implies that the composite operator BF(u) is completely continuous 
from E into E. 

The estimates of Lemma 17 imply that all solutions of equation (38) 
satisfy 

(41) II u ll E < 3. 

Thus, we may take r = 4 in the above theorem. Note that u = h is the 

unique solution of u = A(u, 0) and Hull < 4. Thus, we have proved 

E 

Theorem 42 : The boundary value problem (11) - (13) has a solution for 

all R < 0. 

Skala and Wang (1978) have shown that this is the only solution. 

V. NUMERICAL RESULTS 

Almost all numerical results reported in the past three decades on the 
Berman's equation were obtained as follows: 

First, one makes changes of variables 

(43) f (n ) = g(t) , bn = t. 


Then (11) - (13) become 
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(44) g'~ + g' 2 - gg~ = k 0 = 

b 

(45) g(0) = g"(0) = 0, 

(46) g'(b) = 0; bg(b) = k, 

where b and R are treated unknown a priori. By suitable choice of 
g'(0) and g'"(0) together with (45), one can solve the fourth-order 

equation (by differentiating (44) once) step by step until (if possible) 
g'(b) = 0 for some b. Then R is determined by R = bg(b). Even though 
this method has been popular since Terrill (1964), it is not natural in the 
sense that, in practice, we are given R — the Reynolds numbers, and to find 
the velocity fields. Since from the estimates on a, we can expect that we 
must have a very good control on g'(0) and g'"(0) in order to get 

accurate R from Terrill's method. Previous authors did not give sufficient 
information about how large is b for different R. 

We ran the boundary value problem (14), (12), and (13) under NOS. 2.3 at 
the NASA Langley Research Center by using PASVART, a commercial subroutine to 
solve multipoint boundary value problem by collocation and continuation. The 
results of f, f', f" for R between 0 to -8080 are given in Figs. 3-5. 

It is interesting to note that the asymptotic results obtained by Yuan 

and Finkelstein (1956) as R + -°° are in very good agreement with our 

numerical results as can be seen from Figs. 6-8. 

The perturbed solution for large negative R obtained by Yuan and 


Finkelstein (1956) is 
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(47) f(n) = sin i n + I {cos £ n(1.438 + li/ n — ±_) + 

2 R 2 16 1 sin i 

2 

(0.662 + ~ In | tan j n|)(^- sin | n - ncos |n) - — * - ^ 2 ^ sin y n}. 


Thus, 


“ -? + 0 4 >- 

2 i 

f"(l) = - J- + 0(i), as R > - co. 

8 --ir +0 4>- 


ACKNOWLEDGEMENT 

I wish to thank Dr, Hussaini of ICASE for introducing to me the paper of 


Skalak and Wang (1978) 



-15- 


REFERENCES 

Berman, A. S. (1953), "Laminar Flow in Channels with Porous Walls," J. Appl. 
Phys., _24_, pp. 1232-1235, 

Raithby, G. (1971), "Laminar Heat Transfer in the Thermal Entrance Region of 
Circular Tubes and Two-Dimensional Rectangular Ducts with Wall Suction and 
Injection," Int. J. Heat and Mass Transfer, 14 , pp. 223-243. 

Robinson, W. A. (1976), "The Existence of Multiple Solutions for the Laminar 
Flow in Uniformly Porous Channel with Suction at Both Walls," J. Engrg. 
Math., 10 , pp. 23-40. 

Sellars, J. R. (1955), "Laminar Flow in Channels with Porous Walls at High 
Suction Reynolds Number," J. Appl. Phys., 26 , pp. 489-490. 

Skalak, F. M. and Wang, C. Y. (1978), "On the Nonunique Solutions of Laminar 
Flow Through a Porous Tube or Channel," SIAM J. Appl. Math., 34 , pp. 535- 
544. 

Terrill, R. M. (1964), "Laminar Flow in a Uniformly Porous Channel," 
Aeronautical Quart., 15 , pp. 299-310. 

Yuan, S. W. and Finkelstein, A. B. (1956), "Laminar Pipe Flow with Injection 
and Suction Through a Porous Wall," Trans. ASME, 78, pp. 719-724. 



6 - 



Figure 3 







Figure 4 





8 - 



Figure 5 
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